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The classical Heisenberg model is applied in a Monte Carlo study to investigate the distance 
dependence of the indirect nearest neighbor (NN) exchange and next-nearest neighbor (NNN) su- 
perexchange interaction in EuO, EuS, EuSe and EuTe. For this purpose, first, the dependence of 
the magnetic ordering temperature, i.e., Curie, respectively, Neel temperature for ferromagnetic and 
antiferromagnetic ordering on the exchange constants was determined. This was then employed for 
the analysis of experimental data of hydrostatic pressure experiments. It is shown that all experi- 
mental findings, i.e., the strong increase of the critical temperatures, as well as the transition from 
antiferromagnetic to ferromagnetic ordering for EuTe and EuSe with decreasing lattice parameter 
is well described by the magnetic Griineisen law, in which the exchange constants depend on the 
interatomic distances of the Eu ions in the form of a power law. According to these calculations, 
the indirect NN exchange is characterized by a Griineisen exponent of approximately 20 and the 
NNN superexchange by an exponent of about 10 for all four europium monochalcogenides. The 
latter agrees with Bloch's empirical 10/3 law for the volume dependence of superexchange interac- 
tions in insulating magnetic materials. The Monte Carlo calculations also yield significantly revised 
exchange constants for unstrained bulk material because spin fluctuations at non-zero temperatures 
are taken into account. The strong increase of the exchange constants with decreasing lattice pa- 
rameter provides room for increasing the Curie temperatures in strained epitaxial structures, which 
is important for device applications. 



PACS numbers: 75.10.Hk, 75.30.Et, 75.40.Mg, 75.50.Pp 

I. INTRODUCTION 

The europium monochalcogenides (EuJf, with X = O, 
S, Se or Te) are wide band gap magnetic semiconductors 
with cubic rocksalt crystal structure and increasing lat- 
tice constant as X changes from O to Teii 2 -^. They are 
considered to be model substances for Heisenberg mag- 
nets with spin ordering dominated by indirect nearest 
neighbor (NN) exchange Ji and next-nearest neighbor 
(NNN) superexchange J2 acting between the S = 7/2 
localized magnetic moments of the Eu 2+ ions with half 
filled 4/ shells^. Depending on the sign and magni- 
tude of the exchange integrals J\ and J2 , the EuXs ex- 
hibit different magnetic phases below the critical order- 
ing temperature 5 -. EuO^ and Eu££ are ferromagnets and 
EuTe^ is an antiferromagnct. EuSe is at the borderline 
between ferromagnetic and antiferromagnetic ordering. 
Thus, it shows mctamagnetic behavior— i 10 : 11 , which is 
influenced by additional contributions from dipolar in- 
teractions and crystalline field anisotropics. 

The EuA"s show several outstanding properties, which 
makes them an interesting class of materials, both aca- 
demically and for device applications. In external mag- 
netic fields they exhibit a giant spin-splitting of the 
conduction band and, consequently, extraordinary large 
magneto-optical effects. EuSe shows the largest effective 
g-factor— of up to 18 000 and EuTe the largest magnetic 
field induced energy shifts of the intcrband transitions 13 



observed in semiconducting materials. Potential ap- 
plications are spin-filter devices based on Eu O 14 : 15 : 16 , 
FuS^— — — — — or EuSe^ tunnel junctions, which 
provide spin-polarized electrons due to different barrier 
heights for electrons in different spin states. Also, a huge 
Faraday rotation is observed in EuX o 26 : 27 : 28 : 29 due to 
the spin-splitting of the bands, which results in differ- 
ent refractive indices for left and right circular polarized 
light. Therefore, EuS/EuF 2 and EuSe films have been 
used for high resolution magneto-optical imaging of the 
flux distribution in superconductors^ .. Recent work has 
also demonstrated that EuO can be epitaxially grown 
on silico n 31 : 32 and GaN^ 2 ., which opens new possibilities 
for device realization. Since the Curie temperature of the 
EuA"s can be drastically enhanced by dopin g 32 : 33 : 34 , EuO 
might even become a candidate for practical spintronic 
device applications. 

Introducing strain, either omniaxially through hydro- 
static pressure^^^^^ 2 -^ 1 ^^^ or biaxi- 
ally through epitaxial strai n 11 : 49 : 50 : 51 : 52 : 53 leads to dras- 
tic changes in the ordering temperatures in the EuX com- 
pounds and in some cases even to transitions to differ- 
ent kinds of magnetic ordering. For EuO, hydrostatic 
pressure was found to increase the ferromagnetic order- 
ing temperature Tc from 69 to above 200 KM- and for 
EuS from 16 to almost 180 K 4 ^. Metamagnetic EuSe is 
transformed to a stable ferromagnet already at moderate 
hydrostatic pressures above 0.5 GPa^°- and at higher pres- 
sures T c increases from 4.7 to 70 K at 15GPa-£. EuTe 
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remains antiferromagnetic up to 9 GPa with nearly con- 
stant Neel temperature TV ~ 10 K, but then becomes fer- 
romagnetic with a Tc increasing up to 28 K when reach- 
ing 17GPai£. 

The variations in the magnetic properties of the EuA 
compounds are obviously related to the dependence of 
the exchange integrals J\ and Ji on the interatomic dis- 
tances in the crystal lattice. Already by early theoretical 
work, the basic trend of the EuA compounds from an- 
tiferromagnetic (EuTe) to ferromagnetic ordering (EuS 
and EuO) was attributed to a strong increase of the fer- 
romagnetic NN exchange J\ with decreasing lattice con- 
stant from a = 6.598 A for EuTe to 5.144 A for EuO. For 
the latter, the magnetic ordering is thus dominated by 
the positive NN exchange Ji, whereas for antiferromag- 
netic EuTe the negative NNN exchange J2 dominates. 
Application of hydrostatic pressures p up to 20 GPa, pro- 
duces similar changes in the lattice parameter of up to 
8% compared to the normal bulk values. As a result, 
large changes in the ordering temperatures are induced 
as "- C n 35 - 36 - 37 - 38 - 39 - 40 - 41 - 42 - 43 - 44 - 45 - 46 - 47 - 48 

To derive the dependence of the exchange integrals J\ 
and J2 on the interatomic distances from hydrostatic 
pressure experiments, previous works have employed 
the mean field approximation (MFA) for analysi o 46 ! 48 ' 54 . 
Based on the observation that the Neel temperature Tn 
of antiferromagnetic EuTe does not change appreciably 
under applied pressure, it was reasoned that the NNN ex- 
change Ji is constant in all EuA compounds. Therefore, 
the changes in the magnetic properties were attributed 
solely to changes in the NN exchange J\ and, from the 
simple mean field expressions, a distance dependence of 
Ji(a) was deduced from the observed changes of Tc(p). 
However, it is well known that the mean field approxima- 
tion is exact only at zero temperature, i.e., for prediction 
of the ground state of the system. In particular, the 
mean field model vastly overestimates the magnetic or- 
dering temperatures due to neglection of spin fluctuations 
at finite temperatures. As a result, the mean field ap- 
proximation not only predicts false critical exponents at 
the phase transition but also much underrated exchange 
constants. 

In the present work, we have employed the Monte 
Carlo method (MC) to calculate the magnetic phase dia- 
grams and ordering temperatures of the EuA compounds 
as a function of the exchange interactions. The MC 
method takes the mutual interactions between all spins 
into account and allows for spin fluctuations at T above 
zero. Using finite size scaling techniques, MC correctly 
predicts the transition temperatures and the behavior 
of the order parameters at criticality for a given model 
Hamiltonian 5 - 5 -. Moreover, due to the spin fluctuations 
near the transition temperature in the MC calculations, 
the Neel point of antiferromagnetic ordering as in EuTe 
is found to depend on on the exchange constants J\ and 
Jz 5 ^, in contrast to the mean field approximation, where 
the Neel point depends on J2 only. Thus,the basic as- 
sumption of the previous analyses does not hol d 46 i 48 ' 54 . 



The experimental data for EuAs under hydrostatic pres- 
sure is therefore reexamined by Monte Carlo calculations 
and by considering magnetic Gruneisen power laws 5 -! for 
the distance dependence of the NN and NNN exchange 
constants. We show that for the whole family of EuA 
compounds, the distance dependence of the exchange 
interactions J\{r\) and Jiiyz) can be consistently de- 
scribed by unique Gruneisen exponents of n\ ~ 20 and 
ri2 ~ 10 by which the whole set of experimental data 
of the EuA compounds under hydrostatic pressure can 
well be explained. The obtained exponent of ?i2 ~ 10 for 
J2 is also consistent with Bloch's empirical 10/3 law 5 -! 
for the volume dependence of the superexchange interac- 
tion J ~ j/^ 10 / 3 observed for a wide variety of insulating 
magnetic material systems. 

The paper is organized as follows: In Scc.[Tl]we present 
the model Hamiltonian and briefly discuss the technical 
details of the Monte Carlo calculations. In Sec. IIII1 the 
method is applied to bulk EuTe under ambient pressure, 
demonstrating that the whole magnetic phase diagram 
can be well described and that the exchange constants J\ 
and J2 obtained from Monte Carlo calculations strongly 
differ from literature values derived by the mean field 
approximation. In Sec. IIV1 the method for determina- 
tion of the distance dependence of the exchange con- 
stants is described and applied to EuTe. Due to the pres- 
sure induced transition between antiferromagnetic and 
ferromagnetic ordering, accurate dependencies for both 
Ji(ri) and ^(72) are determined. The approach is then 
extended to EuO, EuS and EuSe in Sees. [V] and IVTl 
revealing that the same functional behavior, i.e., the 
same Gruneisen exponents provide an excellent descrip- 
tion of the experimental data for all EuA compounds. 
In Sec. IVIII the results are compared in detail and the 
applicability of other types of functional dependence of 
Ji(Yi) and ^(^2) discussed. 

II. DETAILS OF THE CALCULATION 

For the calculation of the magnetic properties of the 
EuAs we employed the classical Heisenberg model with 
nearest and next nearest neighbor exchange interaction 
taken into account. The corresponding model Hamilto- 
nian reads as 

H = - J~] JjjSjSj - gfiBH y~] Si , (1) 

where H denotes the external magnetic field and J^- + 
Jji is the total exchange interaction between two spins 
located at lattice sites i and j and 

( Ji : i is NN of j 
Jij = I J 2 : iis NNN of j . (2) 
[ : else 

In the Monte Carlo calculations we considered rhom- 
bohcdral fee clusters of classical spins, where all clus- 
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tcr boundaries are (111) lattice planes. This is a con- 
venient choice of geometry, since antiferromagnetic and 
ferrimagnctic ordering in EuSe and EuTe is comprised of 
ferromagnetic (111) planes and epitaxial EuTe and EuSe 
samples are usually grown in (111) orientatio n 11 ' 12 ' 13 ' 58 . 
The choice of geometry has, however, no influence on the 
results of our calculations. Clusters of up to 32 3 spins 
with periodic boundary conditions were considered. Dur- 
ing a single Monte Carlo step, random orientations are 
generated for every single spin, which are then accepted 
or rejected according to the Metropolis criterion 5 -^. Ob- 
serveables like the total energy, the overall magnetiza- 
tion, the magnetization in the direction of the external 
field, the transverse magnetization and the corresponding 
staggered magnetic moments are computed after every 
Monte Carlo step. Simulations were performed with up 
to N = 10 5 iterations and additional Monte Carlo steps 
for equilibration at the beginning of every run. For a sin- 
gle simulation the temperature T, the external magnetic 
field H and the number of spins are constant. 

The expectation value for the total energy is given by 



Sill 7-Li exp (-0Hj 
■pHi) 



(H) 



Ei=i exp(- 



(3) 



where (3 = l/ks T. The Metropolis algorithm causes the 
total energy to be distributed according to Boltzmann's 
law. Hence, the expectation values for the total energy 
and for the magnetic moment become arithmetic mean 
values in the simulation 



(H) 



1 N 



; = 1 



A' 



(4) 



(5) 



next-nearest neighbor exchange ( J2 = 0) between classi- 
cal spins in the fee lattice, we obtained a critical temper- 
ature T c defined by 2 Ji S 2 /(k B T c ) = 0.3149 ± 0.0008 
in 2 x 10 5 Monte Carlo steps. The prefactor of 2 stems 
from the fact, that according to the definition of Eqn. ([I]) 
the exchange interaction between pairs of spins is always 
added twice. Despite the simplicity of our approach com- 
pared to other sophisticated Monte Carlo routine s 61 ' 62 , 
our result is in excellent agreement the theoretically 
predicted and generally accepted value of Ritchie and 
Fisher^ , of 0.3147 ± 0.0001 for this relation. 

If the NNN exchange interaction is antiferromagnetic 
and IJ2I > Ji, the simulation generates a classical Neel 
state with eight ferromagnctically ordered sublattices as 
the ground state 6 - 3 -. This is a consequence of the antiferro- 
magnetic ordering degenerating into the four equivalent 
(111) directions 5 ^ in the simulation. In this case, four 
pairs of antiferromagnetically aligned sublattices can ro- 
tate freely and the fourth order cumulants are defined 

Q£1 63 
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for zero field and 
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(7) 



(8) 



for nonzero external field. Here, M st and MJ* denote the 
staggered (transverse) magnetization, with transverse re- 
ferring to the component of the magnetization perpendic- 
ular to the external field. 



M st = ^ |M ( 



(9) 



To determine the critical ordering temperatures the 
fourth-order cumulant of the corresponding order param- 
eter is used, which for ferromagnetic ordering at zero ex- 
ternal field is defined as 5 - 5 - 



17 = 1- 



(M 4 



3(M 2 



(6) 



where (M 2 ) and (M 4 ) denote the second and fourth order 
moments of the probability distribution of the magneti- 
zation. The fourth order cumulants show universal values 
at the critical temperature. Thus, during simulation U 
is generated as a function of temperature and recorded 
for various cluster sizes. The different U(T) curves cross 
in a single point at the critical temperature, as shown in 
detail in Sec. IIIIl The transition temperatures obtained 
from our calculations were compared to theoretical pre- 
dictions from high-tcmperature series expansions 6 - -. For 
ferromagnetic nearest neighbor exchange ( J\ > 0) and no 



Mf 



E 

i=l 



M 



(10) 



Eqns. ([9]) and ([10]) are sums of the absolute values of the 
(transverse) magnetization over the eight possible sublat- 
tices. Due to additional anisotropies and/or dipolar cou- 
plings, this degenerate AFM state is not observed experi- 
mentally, but only domains with completely ferromagnet- 
ically ordered Eu (111) planes, which is one possible case 
in our MC simulations. Here it should be also noted that 
in a quantum mechanical treatment it has been showr*^ 4 - 
that a classical Neel state is not an eigenstate of the sys- 
tem. However, Anderson 6 - 5 - showed that the upper limit 
of the error introduced by utilizing the classical Heisen- 
berg model is 1/(Z S), where Z is the number of nearest 
neighbors. Since the Eu-ions carry a relatively large spin 
of S = 7/2, the error in the ground state energy is smaller 
than 2.4%. Thus, a classical treatment considering con- 
tinuously rotating spin vectors is well justified. 



4 



III. EXCHANGE INTERACTIONS IN EuTe AT 
AMBIENT PRESSURE 

In our Monte Carlo study of the exchange interactions, 
EuTe is chosen as test material. This is because de- 
tailed experimental data is available for the H-T phase 
diagram^, which allows a direct determination of the NN 
and NNN exchange interaction J\ and 3i based on the 
antiferromagnetic ordering temperature Tn and the criti- 
cal held at zero temperature Hc(T = 0). In addition, the 
ferromagnetic (Ji) and the antiferromagnetic exchange 
(J2) are relatively balanced in EuTe. Therefore, pro- 
nounced changes in magnetic ordering occur when hy- 
drostatic pressure is applied. 



A. Experimental Results 

The magnetic properties of EuTe were determined 
using DC and AC SQUID magnetometry measure- 
ments of high quality 4^tm thick (lll)-oricntcd epitax- 
ial layers grown by molecular beam epitaxy on BaF2 
substrate s 13 ' 66 ' 67 . As a result, the magnetization and AC 
susceptibility was obtained as a function of both temper- 
ature and external magnetic field up to 7T. The exter- 
nal field direction was applied in the (111) growth plane, 
which is also the easy plane of the magnetization. 

Experimental magnetization curves M(T) at different 
applied external fields from 5.5 to 7 T are shown in 
Fig. [lja). The magnetization curves exhibit clear peaks 
at the phase transition between the antiferromagnetic 
and the paramagnetic phase. High magnetic fields were 
applied in order to obtain information about the critical 
field at zero temperature Hc(T = 0). Field-dependent 
AC susceptibility curves were also measured. As shown 
in Fig. [ljb), for T = 5K the susceptibility is essentially 
constant below the critical field, corresponding to a lin- 
ear increase in the magnetization. This arises from the 
continuous reduction of the relative angle between the 
spins in the adjacent (111) lattice planes from 180° in 
the AFM II state at zero field to almost zero at the crit- 
ical field He, as illustrated schematically by the arrows 
in Fig. [He). Thereby the spin orientation changes from 
initially perpendicular to the external field to finally par- 
allel to the external field direction when H reaches the 
critical field. The critical field is thus given by the dis- 
continuous drop in the susceptibility at this point (see 

Fig. up))). 

Figure [He) displays the susceptibility as a function of 
temperature at zero external field, showing a broad peak 
with a maximum slightly above 10 K. However, the phase 
transition from antiferromagnetic to paramagnetic corre- 
sponds to the peak in the specific heat, which according 
to Fisher's rule^ coincides with the maximum slope, i.e., 
the inclination point in the AC susceptibility. As indi- 
cated by the arrow in Fig. [He) , the phase transition thus 
occurs at T N = 9.85 ± 0.05 K. 

Figure^d) shows the phase diagram of EuTe compiled 
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FIG. 1: Magnetic properties of bulk EuTe at normal pressure 
derived from SQUID measurements and Monte Carlo simula- 
tions, (a) Temperature-dependent magnetization M(T) mea- 
sured for various in plane external magnetic field values from 
5.5 to 7 T. (b) AC susceptibility vs. external magnetic field at 
5 K indicating the critical field as discontinuity at B « 6.2 T. 
(c) AC susceptibility \ measured as a function of temperature 
at zero external field. The inclination point indicated by the 
arrow yields a Neel temperature of 9.85±0.05 K. (d) Deter- 
mination of the Neel temperature from Monte Carlo data us- 
ing the fourth order cumulant of the staggered magnetization 
U(M st ) for different system sizes. The curves show a com- 
mon intersection at the Neel temperature, (e) Phase diagram 
of EuTe: Symbols with error bars correspond to maxima or 
inflection points in M(T) for high, respectively, small exter- 
nal fields; squares to the discontinuity in x(^) or the inflec- 
tion point in x(T)- The solid line indicates the experimental 
Hc(T) phase boundary and the dash-dotted line represents 
the boundary obtained by the Monte Carlo calculations. The 
T 3//2 extrapolation (dashed line) of the experimental data to- 
wards T — yields a critical field of He = 7.56 T. 



from the experimental data (symbols and solid line), 
where the thick solid line indicates the experimentally 
determined Hc(T) phase boundary. A T 3 / 2 curve, which 
according to spin wave theory is the low-temperature be- 
havior of the critical field, was fitted to the experimental 
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Reference 


Exper. results / Analysis 


Ji (K) 


J 2 (K) 


Oliveira^ 


H c , 1 MFA 


0.100 


-0.215 


Zinni 


H c , 6 1 MFA 


0.060 


-0.200 


Wachtes^ 


H C ,T N 1 MFA 


0.043 


-0.150 


Kobles^ 


He, Tn 1 MFA 


0.060 


-0.160 


Kunes^ 


ab initio LDA+f/ (U = 6eV) 


0.110 


-0.320 


Our study 


H C ,T N / MFA 


0.035 


-0.156 




H c , T N 1 MC 


0.192 


-0.313 



TABLE I: Comparison of the exchange constants Ji/fcs and 
JiI^b of EuTe determined by the analysis of experimen- 
tal data for the Neel temperature Tn, the critical field He 
or the paramagnetic Curie temperature 9c using the mean 
field approximation (MFA) or the Monte Carlo (MC) method 
(present work). Also listed are the exchange constants derived 
by Kunes et al. from ab initio calculations. 



critical points between 2K and 6K, yielding 
H C (T) = H c (0) [l - A(T/XV) 3 / 2 



(11) 



Extrapolation of the the measured Hc(T) to T = thus 
yields the critical field at zero temperature of -ffc(O) = 
7.56 ± 0.02 K. The dashed line in Fig. fjje) represents 
Eqn. dTTJ) with the coefficient A = 0.50±0.01. The exper- 
imental data as well as TV, Hc(0) and A are in excellent 
agreement with previous results of Oliveira et al£. 



B. Exchange constants and phase diagram from 
MC calculations 



In most previous studies^&Si, the mean field anal- 
ysis was used to determine the exchange integrals in 
EuTe, because it provides simple analytic expressions for 
the critical field at zero temperature Hc(0), the Neel 
temperature TV as well as the paramagnetic Curie tem- 
perature 6c as a function of Ji and Ji. For type II- 
antifcrromagnctic ordering in an fee lattice with NN and 
NNN exchange interactions, the mean field approxima- 
tion (MFA) yields 



itMFA 



(0) 



MFA 
N 



-4S(6Ji + 6J 2 )/(<?/i B ) 
^S(S + l)(-6J 2 )/k B 



e 



MFA 



-S(5 + l)(12Ji + 6J 2 )/* 



(12) 
and (13) 

. (14) 



where g = 2 and S = 7/2 for the magnetic moment of the 
Eu 2+ ions. Inserting our experimental values for Hc(0) 
and Tn into Eqns. (fT2")) and p3[) and solving for Ji and J2 
yields J^ FA /k B = 0.035 K and Jf FA /k B = -0.156 K. As 
shown in Tab. [H the values are consistent with previous 
mean field studies, especially those, which applied the 
same analysis of the experimental Hc{0) and Tjv based 
on MFA Eqns. ((I2J| and $F$. 




-0.1 -0.15 -0.2 -0.25 -0.3 -0.35 -0.4 
Next-nearest neighbor exchange J 2 / k B (K) 

FIG. 2: Neel temperature of EuTe as a function of the 
nearest neighbor exchange constant Ji under the constraint 
(Ji + J2)/fcs = — 0.12 IK, which is a consequence of eval- 
uating the experimental critical field at zero temperature 
He(T = 0) — see in the text for details. We obtain a linear 
function (squares and solid line) which is shown in comparison 
with the corresponding mean field relation (dashed line). At 
J 2 = -0.313 K (Ji = 0.192 K) the Monte Carlo curve reaches 
the experimental Neel point of 9.85 K. 



In the Monte Carlo calculations, the transition tem- 
perature Tjv(0) (Tn(H ^ 0)) is deduced from the tem- 
perature dependence of the 4 th order cumulants of the 
staggered (transverse) magnetization U*±y As described 
in Sec. UH these cumulants show universal values at the 
critical temperature TV independent of cluster size. This 
is demonstrated in Fig.fTJJd) for the case of H = 0. Using 
the exchange constants jf IFA and J2 IFA derived from 
the mean field analysis, the Monte Carlo calculations 
yield a Neel temperature of 5.45 K at zero external field, 
which is in strong disagreement with the experimental 
value of 9.85 K. This clearly demonstrates that the ne- 
glection of spin fluctuations in MFA leads to a vast under- 
estimation of the exchange constants, an effect that has 
been already noted in previous theoretical studie s 56 ' 63 . 
In the MC calculations, moreover, the critical Neel tem- 
perature TV is found to depended significantly not only 
on the antifcrromagnctic exchange constant J2 but also 
on the ferromagnetic NN exchange constant Ji, in con- 
trast to the MFA approximation where Tn depends only 
on the antiferromagnetic exchange — see Eqn. (|13[) . This 
is due to the fact that in the mean field approximation, 
for type II antiferromagnetic ordering the Ji exchange 
between the 6 NN Eu 2+ ions within the ferromagnetic 
(111) planes exactly cancels with the J\ exchange to the 
6 NN Eu 2+ ions within the antiferromagnetically coupled 
neighboring (111) lattice planes. This does not apply for 
the MC calculations because of the non perfect antifer- 
romagnetic spin alignment at nonzero temperatures that 
results from spin fluctuations. 

In order to determine the exchange constants from the 
Monte Carlo calculations, we have systematically calcu- 
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lated the Neel temperature as a function of both exchange 
constants J\ and J 2 . As a boundary condition, we take 
advantage of the fact that in the limit of T — 0, the crit- 
ical field He of the Monte Carlo calculations converges 
exactly to the mean field value of Eqn. (|T2|) . This is be- 
cause spin fluctuations are absent at T = 0, and thus 
the MFA represents the exact ground state of the sys- 
tem. Therefore, the experimental value of Hc(T = 0) 
and Eqn. (| 1 2|) determine the sum of the exchange con- 
stants as 

(Ji + J 3 )/fc B = -0.121 K. (15) 

in EuTe. This eliminates one independent variable, i.e., 
with this condition, only J 2 must be varied for the cal- 
culation of Tn. The resulting dependence is plotted in 
Fig. [2] (squares and solid line) . Evidently, T/v varies al- 
most perfectly linearly and can be represented by the 
relation 

^ C | (Ji+j2)/fcB= „ . 121K = 1.22-27.57 J 2 /k B (16) 

within the range —0.35 < J 2 /fcs < —0.15. Solving 
Eqn. p6p for J 2 and inserting the experimental Neel 
point of bulk EuTe Tjy = 9.85 K (horizontal dash-dotted 
line in Fig. [2) yields Jf c /k B = -0.313 K and hence 
J^ c /k B = 0.192 K from Eqn. (15]) as the intrinsic ex- 
change constants of bulk EuTe. It is noted, that calcu- 
lating Tm without the constraint of Eqn. (fTS)) yields a 
function T/v(Ji, J 2 ) that depends nonlinearly on J\ and 
J 2 , in contrast to the MFA Eqn. (fT3"|) . which predicts 
only a linear dependence on J 2 — see Sec. IIVI for further 
details. From the MC calculations, the Neel temperature 
as a function of J\ and J 2 is found to be well described 
by 

T^ c » (-15.3 Ji -40.8 J 2 )/k B 

= (S + 1) (-1.46 Ji- 3.89 J 2 )/k B (17) 

in a linear approximation in the vicinity of the intrinsic 
EuTe exchange constants, demonstrating that the Neel 
temperature indeed depends strongly on both exchange 
constants. 

As demonstrated by Tab. HI which compares our de- 
rived set of exchange constants with previously published 
ones, our values are nearly twice as large as those derived 
from mean field analysis. Thus, by ncglection of spin 
fluctuations the exchange parameters are vastly under- 
estimated. Remarkably, the exchange constants derived 
from our Monte Carlo calculations are in good agree- 
ment with recent ab initio calculations of Kunes et ali^ 
using the local-density-approximation method including 
strong Coulomb repulsion within the 4/ shells (LDA+J7). 
In particular, our NNN exchange constant J 2 , which in 
Ref. [lO depends very weakly on the Coulomb parameter 
U matches the ab initio result very well. 

With the new exchange parameters, we can now calcu- 
late the whole magnetic phase diagram of EuTe using the 
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FIG. 3: (color online) Critical temperatures (left scale) and 
exchange interactions Ji and J2 (right scale) in EuTe plot- 
ted as a function of the lattice constant a. Open and closed 
symbols (circles^, squares^, triangles^), correspond to Neel 
(Tat), respectively Curie points (7b) determined from hydro- 
static pressure experiments. Solid lines: Tjv(a) and Tc(a) de- 
rived from Monte Carlo calculations using a Ji(a) and J2(a) 
dependence given by the Griineisen law (|18[) with best fitting 
power law exponents of n\ — 20.6 ± 0.4 and ri2 = 10.4 ± 0.5, 
respectively. The corresponding variation of Ji(a) and J2(a) 
as functions of the lattice parameter is shown as the dashed 
and dash-dotted line, respectively. 



4 th order cumulants of the staggered (transverse) magne- 
tization for various external magnetic fields and cluster 
sizes. The resulting phase boundary Tff(H) is depicted 
as dash-dotted line in Fig. [TJc). As expected, the calcu- 
lated Hc(T) approaches the experimental value of 7.56 T 
in the limit of T — > and the calculations nearly follow 
the experimental Hc{T) boundary. The fact, that at 
T > the calculated Hc(T) values are slightly lower 
than the measured ones and that at low temperatures, 
the calculated critical field varies linearly with tempera- 
ture instead of obeying a T 3 / 2 behavior is a well known 
consequence of applying a classical (S = 00) instead of 
the quantum mechanical S = 7/2 model in our calcu- 
lations. As already noted in Sec. HH the error in the 
ground state energy introduced by this simplification is 
of the order of less than 2.4 % for our type of system. 



IV. DISTANCE DEPENDENCE OF EXCHANGE 
INTERACTIONS IN EuTe 

Under hydrostatic pressure, the EuTe lattice constant 
can be compressed from its normal bulk value of an = 
6.589 A to about 6.15 A at a pressure reaching 17GPa4*. 
This corresponds to a 7% reduction of the lattice con- 
stant and of the interatomic distances of the Eu 2+ ions in 
the crystal , where in the fee lattice of EuTe the NN Eu 2+ 
distance is r\ = a/\/2 and the NNN distance r 2 = a. The 
resulting changes in the magnetic ordering temperatures 
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obtained by hydrostatic pressure experiment s 38 ! ' are 
compiled in Fig. [31 where the open symbols represent the 
measured Tjv and the full symbols the Tc values plotted 
as a function of lattice constant. Since at ambient pres- 
sure, the antiferromagnetic exchange J 2 of EuTe is larger 
than the ferromagnetic exchange J\ , a type II antiferro- 
magnetic ordering^ occurs below the Neel point of Tjv 
= 9.85 K. With increasing pressure, i.e., decreasing lat- 
tice constant, the Neel temperature remains practically 
constant at Tjv ~ 10 K, but at ss 9GPa or 5% compres- 
sive strain EuTe becomes ferromagnetio 46 ' 47 with rapidly 
increasing Curie Tc that rises up to 28 K at 17GPa 4 ^. 
The observed phase transition from antiferromagnetism 
to ferromagnetism at a = 6.29 A implies that at smaller 
atom distances the NN J\ becomes the dominating ex- 
change mechanism. 

The influence of the inter-atomic distances on ex- 
change constants has been a subject of many theoreti- 
cal studies 1 "^ ' 54 ! 71 ' 72 ! 73 . However, indirect and superex- 
change mechanisms involve complex integrals such that 
up to now no general analytic expressions for their dis- 
tance dependence have been derived theoretically. An 
empirical power law dependence, referred to as the mag- 
netic Gruneisen law, has been proposed by Bloch^ 7 -, i.e., 

J(r) = Jo , (18) 

where Jo = J{ro) and tq are the exchange interaction 
and interatomic distance at normal pressure, and n is 
the scaling exponent. As shown in Ref. l57l . this depen- 
dence well describes the observations for many magnetic 
semiconductors or insulators such as the Mn and Gd 
chalcogcnides or iron oxides, for which the power law 
exponent n shows a universal value of around 10 for the 
magnetic superexchange^ 7 -. This also yields the empirical 
10/3 law for the volume dependence of superexchang 
of J(V) = MV/V0)- 10 / 3 . 

To test if the Gruneisen dependence of Eqn. (fT8"|) ade- 
quately describes the atomic distance dependence of the 
exchange integrals in EuTe, we have performed a series 
of Monte Carlo calculations of the Neel and Curie tem- 
perature as a function of the exchange integrals in order 
to fit the experimental TV (a) and Tc(a) data of Fig. [3] 
using the power law exponents n\ and n 2 for the NN 
and NNN exchange interactions J\ and J 2 as free pa- 
rameters. In these calculations, the ferromagnetic NN 
exchange J\ > and the antiferromagnetic NNN ex- 
change J 2 < were varied independently in the range 
of 0.190 < Ji < 0.73 and 0.315 < -J 2 < 0.615 and the 
corresponding critical temperatures were derived as de- 
scribed in detail in Sees. ITT1 and IIIII For all |Ji| < [ ,7a | , 
the MC calculations yield AFM II ordering, whereas FM 
ordering results for all | Ji | > | J 2 1 . Figures Ufa) and (b) 
show the calculated Tjv and Tc values (open symbols) 
as a function of J\ and J2, respectively. Since we are 
interested in the hydrostatic pressure effect, the small- 
est values of J\ and J2 were chosen close to the exchange 
parameters of bulk EuTe at ambient pressure (filled sym- 



bols) and the maximum values correspond to hydrostatic 
pressures up to about 17GPa. 

As shown by Fig.UJ the Monte Carlo calculations yield 
a strongly nonlinear dependence of Tjv on J\. Hence, 
the calculated data points were approximated by second 
degree polynomials for Tjv(Ji, J 2 ) and Tc(Ji, J2) of the 
general form 

T{J 1 ,J 2 )=A+BJ 1 +CJ 2 +Djf+EJ 1 J 2 +FJ 2 2 . (19) 

These approximations fit the Monte Carlo data with bet- 
ter than ±0.1 K accuracy and are represented as solid 
lines in Figs.|3Ja) and (b). 

Using the normal pressure exchange constants J±.q = 
0.192 K and J 2 , = -0.313 K determined in Sec. UHl the 
two branches of experimental data sets for the antiferro- 
magnetic Tjv(a) and ferromagnetic Tc(a) of Fig. [3] were 
fitted using Eqns. (fT5)) and (|19[) with common exponents 
n\ and n 2 and r\ = a/y(2) and r 2 = a. All experimen- 
tal data points are weighted equally in the least square 
fit routine, which was performed on a logarithmic scale 
since equal weights may cause one branch to dominate if 
there is a difference in the magnitude of the dependent 
variable. 

From this modeling, a Gruneisen exponent of n\ = 
20.6 ± 0.4 for the NN exchange J x and of n 2 = 10.4 ± 0.5 
for the NNN exchange J 2 were obtained. As demon- 
strated by the solid lines in Fig. [3l with these parameters 
the whole body of experimental findings, including the 
approximately constant Neel temperature Tjv at small 
hydrostatic strain, the transition from antiferromagnetic 
to ferromagnetic ordering at a — 6.29 A and the steep 
superlinear increase of the Curie temperature at high 
pressures, i.e., small lattice constants, are exactly repro- 
duced. Moreover, the obtained power law exponent n 2 
for J 2 (r 2 ) is in excellent agreement with Bloch's^ 10/3 
law for the volume dependence of superexchange. 

The such obtained dependence of the NN and NNN 
exchange integrals J\{a) and J 2 (a) as a function of lat- 
tice parameter a is presented in Fig. [3] as dashed and 
dash-dotted line, respectively. Evidently, both exchange 
constants strongly increase with decreasing lattice con- 
stant. However, J\(a) increases much more rapidly than 
J 2 (a) due to the two times larger power law exponent. 
Therefore, the two curves intersect at a = 6.29 A, where 
Ji/ks = —J2/ks = 0.51 K, and at smaller a, the ferro- 
magnetic Ji becomes the dominating exchange mecha- 
nism. For fee lattices with competing ferromagnetic NN 
exchange and antiferromagnetic NNN exchange interac- 
tions, this is exactly the condition for the material to 
become ferromagnetic^. 

Our results are in severe contrast to the previous mean 
field analysis of Goncharenko et al4&, who concluded 
from the negligible variation of Tjv in the antiferromag- 
netic state with changing lattice constant that the NNN 
exchange J 2 in EuA should not depend on the lattice 
parameter. Consequently, the whole variation of Tc(a) 
was attributed solely to changes in J± (a) using the mean 
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FIG. 4: Monte Carlo calculated critical temperatures of a bulk fee Heisenberg system as functions of the strength of (a) the 
nearest neighbor exchange interaction Ji and (b) the next-nearest neighbor exchange interaction J2, showing antiferromagnetic 
to paramagnetic transitions (squares) when Ji < IJ2I and ferromagnetic to paramagnetic transitions (circles) when Ji > \Ji\. 
Solid lines represent polynomial fit functions, which are quadratic in J\ and J2. 



field expression for the Curie temperature of 

Tc FA = \S{S+ 1) (12 Ji + 6 J 2 )/k B . (20) 

with constant Ji for data analysis. On the contrary, our 
calculation show that the broad plateau of Tjv(a) for lat- 
tice constants around a = 6.42 A just results from the 
fact that in the antiferromagnetic phase the ferromag- 
netic exchange drops faster than the antiferromagnetic 
exchange as the lattice constant increases. 



V. EXCHANGE INTERACTIONS IN EuO AND 

EuS 

As shown in the previous sections, the exchange con- 
stants obtained by Monte Carlo calculations strongly dif- 
fer from previously published values. Therefore, to evalu- 
ate the distance dependence of the exchange constants of 
EuO and EuS, first the bulk values under ambient pres- 
sure have to be reexamined by the Monte Carlo method. 



A. Exchange constants at ambient pressure 

EuO and EuS are low temperature fcrromagncts with 
Curie temperatures Tc of 69.15 K-Zi and 16.6 K 74 i 75 , re- 
spectively. The magnetic properties arc determined 
mainly by the dominant ferromagnetic NN exchange in- 
teraction Ji in both materials. Compared to the case 
of antiferromagnetic EuTe, where only the Neel temper- 
ature and the critical field at T = OK are needed to 
deduce the exchange constants, the determination of 3\ 
and Ji for ferromagnetic EuO and EuS is much more in- 
volved. As a result, there exists a substantial variation in 
the reported exchange constants for bulk EuO and EuS 



deduced from different experimental techniques like in- 
clastic neutron scatterin g 74 i 75 i 76 , specific hea t 77 ' 78 i 79 ' 80 , 
nuclear magnetic resonanc o 81 i 82 ' 83 (NMR) and spin wave 
resonance measurements^ (see, e.g., Passell et al^ for 
a review). Especially the values for the NNN exchange 
constant Ji in EuO and EuS differ by up to a factor 
of two in literature and whether Ji is ferromagnetic or 
antiferromagnetic in EuO is still a matter of debate. 
The most recent results based on inelastic neutron scat- 
tering studies on sing le crystals of EuO 76 and EuS 7 ^ 
yielded (J^°/k B , Jf$°/k B ) = (0.625 K, 0.125 K) and 
{Jffl k B, Jif/k B ) = (0.221 K, -0.100 K), respectively, 
consistent with Passell et al.'sZ 4 - analysis on powdered 
samples. Notably, a ferromagnetic NNN exchange inter- 
action was obtained for EuO. 

In all studies, the sum of J = J\ + J 2 has been more 
reliably determined than the individual NN and NNN 
exchange interactions, and this sum is quite consistent 
among the various studies. For EuO single crystals, 
J Eu °/k B = 0.755 K was obtained by Comment et alM 
from NMR measurements, in agreement with neutron 
scattering studies by Mook at ali 7 ^, and this value also 
agrees with the results obtained by neutron scattering 7 - 4 - 
and specific heat measurements^ on powdered samples. 
A very good agreement for J± + Ji also exist among re- 
spective studies for Eu S 74 ' 75 i 80 ' 81 , from which we calcu- 
late J//c| uS = 0.121 ± 0.003 K as mean value. 

To determine the exchange constants by the Monte 
Carlo method, we again performed a series of calcula- 
tions for the model Heisenberg Hamiltonian of Eqn. (Q]) 
with the NN and NNN exchange interactions J\ and Ji 
varied independently over a wide range of 0.5 < J\ < 2.5 
and 0.4 < — Ji < 0.8. We find ferromagnetic ordering for 
all combinations of Ji > — Ji and determined the corre- 
sponding critical ordering temperature as a function of 
J\ and Ji as presented in Fig. [5] (a) and (b), respectively. 
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Evidently, Tc increases linearly with increasing NN ex- 
change constant Ji, but decrease when J 2 increases. As 
shown by the solid lines in Fig. [5j the general dependence 
of Tc on the exchange constants can be well described 
by the relation 

T^ c = 79.0 Ji/k B + 55.9 J 2 /k B (21) 

= 0.627 ^-S (5 + 1) 12 (Ji + 0.708 J 2 )/k B , 
o 

which differs considerably compared to the mean field 
expression of Eqn. (|2"0|) . Yet, the critical coupling 
J\/{k B Tc) closely resembles the theoretical predictions 
from high-temperature series expansion^. Inserting the 
EuO and EuS exchange values of MookZ& and Bohn et 
ali^ in our Monte Carlo relation of Eqn. (f2"T]) yields Curie 
temperatures of only 56 K and 12 K, respectively, which is 
much lower than the measured experimental values. This 
shows that, like for EuTe, the exchange constants have 
been considerably underestimated in both materials. 

Using the experimental values for J = J\ + J 2 quoted 
above, the critical Curie temperature Tc can be calcu- 
lated MS (X function of the exchange constant J\ using 
Eqn. (|2"Tj) and J 2 = J — J\. The results arc plotted as 
solid lines in Fig.^a) and (b) for EuO (J Eu ° = 0.755 K) 
and EuS (J EuS = 0.121 K), respectively. From the in- 
tersection of these lines with the respective experimen- 
tal T c values of 69.15 KZi and 16.6K2iiS (horizontal 
dashed lines in Fig. O, the bulk exchange constants of 
J Eu ° = 1.169K and J 2 Eu ° = -0.414K are obtained 
for EuO and of J EuS = 0.427 K and J 2 EuS = -0.306 K 
for EuS. In Fig. [5J also plotted are the Tc values ex- 
pected from the mean field approximation (dash-dotted 
line, Eqn. (HOI)) as well as from a series expansion es- 
timate proposed by Passell et al.— (dashed line) given 
by 

T SE, cst . = o.790 ^ 5 (S+l) 12 (Ji + 0.619 J 2 )/A:b ■ 

(22) 

Evidently, in both cases, much higher critical tempera- 
tures are predicted for a given set of exchange constants, 
i.e., from the observed transition temperatures, the ex- 
change integrals are strongly underestimated. Moreover, 
a ferromagnetic NNN exchange would be suggested for 
EuO. Our results exclude such a ferromagnetic NNN ex- 
change, i.e., a negative J 2 is obtained for EuO even if 
jEuO j g underestimated by as much as 15 %. Thus, as 
already found for the EuTe case, our Monte Carlo calcu- 
lations greatly revise the bulk exchange constants. 



B. Distance dependence of EuO and EuS exchange 
constants 

For EuO and EuS, the Curie temperature Tc strongly 
increases with increasing hydrostatic pressure, i.e., de- 
creasing lattice constant. This is illustrated by Figs. [7| 
and [SI where the experimentally determined Tc values of 




Nearest neighbor exchange constant Jj (K) 



(b) 


200 










g 




□ 


B 


J i =2.5K 




£ 


150 






H B . 


■ □ 


3 




□ 


□ 






C3 








Q □ 


□ 


<D 
ft 


100 










E 
B 












■£ 


50 


□ 


□ 


D □ 




u 





□ 


J,= 


0.5 K 


□ 
















-0.4 


-0.5 


-0.6 -0.7 


-0.8 



Next nearest neighbor exchange constant J 2 (K) 

FIG. 5: Monte Carlo calculated critical temperatures of a 
bulk fee Heisenberg system as functions of the strength of 
(a) the nearest and (b) next-nearest neighbor exchange inter- 
actions (Ji and J2, respectively), showing ferromagnetic to 
paramagnetic transitions if Ji > | ^"2 1 . The parameter range 
covers EuO and EuS under hydrostatic pressures between 
and 20 GPa. Solid lines represent fit functions, which are lin- 
ear in Ji and J 2 — see Eqn. (|21[) . 

Eu0^4^ and EuS 48 are pi otte d as a f unc tion of the lat- 
tice constant. At hydrostatic pressures around 20 GPa, 
corresponding to a 6-8% reduction of the lattice con- 
stant, T c is as high as 200 and 180 K^, respectively. 
For EuO, the experimental data Tc(a) of Fig. [7] is com- 
piled from three independent investigation a 35 ' 44 i 45 , with 
McWhan et al.'s data combined with the pressure- volume 
relation taken from Ref. [8^. For EuS, the data is taken 
from Ref. [4SL For EuO, experiments at even higher hy- 
drostatic pressures up to 31 GPa^ have revealed that the 
ferromagnetic ground state becomes unstable at around 
23 GPa and that the Curie temperature drops sharply 
afterwards instead of further increasing. In Ref. |45| this 
behavior was attributed to sf hybridization competing 
with sf exchange in this pressure range. Therefore, we 
restrict our analysis to the 0-20 GPa range, i.e., lattice 
constants above 4.9 A, where such effects seem not to be 
of importance. 

To determine the interatomic distance dependence of 
the exchange constants, we proceed in the same man- 
ner as described in Sec. IIVI by fitting the calculated 
T,¥ c (Ji, J 2 ) dependence of Eqn. obtained by the 
Monte Carlo calculations to the data of the hydrostatic 
pressure experiments, applying the magnetic Griineisen 
law (Eqn. I|18p) as functional dependence for the NN and 
NNN exchange constants. As input parameters we use 
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FIG. 7: (color online) Critical temperatures and exchange in- 
teractions in EuO as functions of the lattice parameter; filled 
symbols (circles^, triangles^, squares^) correspond to ex- 
perimental Curie points; the solid line represents the least 
square fit of Monte Carlo critical temperatures Tc(Ji,J2) 
based on the magnetic Griineisen law Ji(ri) ~ r~ n ' . The cor- 
responding dependence of the nearest (Ji) and next-nearest 
neighbor exchange interaction (J2) on the lattice parameter 
is plotted as the dashed and dash-dotted line, respectively. 
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FIG. 8: (color online) Critical temperatures and exchange 
interactions in EuS as functions of the lattice parameter; 
squares^ correspond to experimental Curie points; the solid 
line represents the least square fit of Monte Carlo critical tem- 
peratures Tc(Ji, J2) based on the magnetic Griineisen law 
Jiifi) ~ r~"\ The corresponding dependence of the nearest 
(Ji) and next-nearest neighbor exchange interaction (J2) on 
the lattice parameter is plotted as the dashed and dash-dotted 
line, respectively. 



the bulk exchange constants J|g and J20 determined 
in the previous section and treat the power law exponents 
n\ and 712 in Eqn. (|18|) as adjustable parameters. It turns 
out, that because the NN exchange 3\ in EuO and EuS is 
always much larger than the NNN exchange J2, the Tcifl) 
dependence is quite insensitive to the variation of J2 as a 
function of lattice constant, i.e., the fit yields only unreli- 
able values for 712 . Because for EuTe we have already con- 
firmed Bloch's 10/3 law for the volume dependence of the 
NNN superexchange integral J2, we have therefore cho- 



sen to fix the distance dependence of Ji(a) proportional 
to r^ 10 for EuO and EuS as well. From the fit, we then 
obtain nf u ° = 19.6 ± 0.4 and nf uS = 22.4 ± 0.3 as the 
Griineisen exponents for the NN exchange interaction. 
The resulting lattice constant dependence of the Curie 
temperatures Tc(a) and exchange interactions J\{a) and 
J 2 (a) are plotted in in Figs. [7] for EuO and [8] for EuS 
as solid, dashed and dash-dotted lines, respectively. Ev- 
idently, an excellent fit with the experimental data is 
obtained over the whole lattice parameter range for both 
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materials. This is also an indication that the choice of 
112 = 10 is a reasonable assumption. It is also noted that 
due to the about a factor of two larger Griineisen expo- 
nent of Ji compared to that of J2, at high hydrostatic 
pressures (small lattice constants), the NN exchange J\ 
is as much as five times larger than the NNN exchange 
J2. Thus, the ferromagnetic NN exchange completely 
dominates the magnetic behavior of both materials. 
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FIG. 10: (color online) Critical temperatures and exchange 
interactions in EuSe as functions of the lattice parameter; 
filled symbols — triangles^ - and squares 4 ^ correspond to ex- 
perimental Curie points. The solid line (Tc) represents the 
least square fit of the Monte Carlo Tc(Ji, J2) to the exper- 
imental Curie points based on the magnetic Griineisen law 
Jiifi) ~ rj ni . The corresponding nearest and next-nearest 
neighbor exchange interaction, Ji and J2 as functions of the 
lattice parameter a, are represented as the dashed and dash- 
dotted lines, respectively. The inset shows the AFM I-FM 
transition occurring in EuSe at hydrostatic pressures around 
0.5 GPa. The solid graph (TW), which is based on an extrap- 
olation of Ji (a) and Ji (a) towards do corresponds to a theo- 
retical AFM II Neel temperature which is energetically com- 
pareable to the experimentally observed AFM I Neel point 
(open symbo l 11 ' 40 ). 



VI. EXCHANGE INTERACTION IN EuSe 

Unlike EuO, EuS and EuTe, which exhibit stable 
magnetic low temperature phases, EuSe is a metam- 
agnet with at least four different known ordered mag- 
netic phases, i.e., two antiferromagnetic phases of type I 
(AFM I) and type II (AFM II), a ferrimagnctic phase 
(FiM) and a ferromagnetic phase (FM). Figure [5] shows 
the H-T phase diagram of unstrained EuSe derived from 
susceptibility measurements on several micrometer thick 
epitaxial layers^, and the corresponding different spin 
configurations are illustrated by the arrows. The phase 
boundaries shown in Fig. [9] are in good agreement with 
earlier publications*^. In particular, the AFM I Neel 
point of T)v,afm 1 = 4.7 K and the critical field of 
H c , 3 (0) = 6.05 T for the transition from AFM II to 
FiM obtained by linear extrapolation of the experimental 
AFM II to FiM phase boundary are in excellent agree- 
ment with Refs. andflol. 

The mctamagnctic behavior of EuSe at ambient pres- 
sure shows similarities to the situation in EuTe at the 
AFM II to FM transition that occurs at a hydrostatic 
pressure of 9 GPa (see Fig. [3]) , where Ji ss | J 2 1 . The 
appearance of an AFM I phase at this pressure that ac- 
companies the AFM II to FM phase transition^ shows 
that at this pressure EuTe is metamagnctic too. Thus, in 
both materials a metamagnetic behavior occurs when J\ 



and IJ2I are approximately equal and cancel each other. 
Then otherwise negligible additional interactions come 
into play. As shown by Fig. 1 101 when applying hydro- 
static pressure of 0.5 GPa^ EuSe becomes a stable fer- 
romagnet and at increasing pressures up to 15 GPa, i.e., 
6% reduction of the lattice constant, the ferromagnetic 
ordering temperature increases from 4.7 to above 70 K^. 
The corresponding experimental Tc (a) data of EuSe is 
displayed as filled symbols in Fig. [TUJ Lechner et aljii 
also showed that the introduction of only little biaxial 
strain in EuSe drastically expands the boundaries of the 
AFM II phase and causes the AFM I phase to disap- 
pear completely. Thus, an AFM II to paramagnetic Neel 
point TV,afm 11 is observed in strained EuSe, which is 
again similar to the situation in EuTe. 

The magnetic phase diagram of EuSe and its metasta- 
bility for already small lattice deformations^ cannot be 
described by isotropic NN and NNN Heisenbcrg exchange 
interactions alone. Especially, the ferrimagnetic and the 
AFM I phases, which show a magnetic structure with a 
periodicity of three, respectively, four atomic layers re- 
quire further distant exchange interactions and/or other 
types of magnetic interactions such as long-range dipolar 
interactions^. Moreover, using MFA and a Hamiltonian, 
which includes the exchange interaction up to the third 
nearest neighbor (J3) and dipolar interactions Fukuma 
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et ali^Il showed that the critical field Hc,3 at T — is in- 
dependent of J 3 and depends only on the sum (Jj + J 2 ) 
as well as the dipolar coupling strength D. A good esti- 



2(Ji l0 + Ja,o) = -3^9Mb^c,3(0) - 



Here D XX (Q L ), D xx (0) and D xx (2Q L /3) correspond 
to the dipole coupling strength for spins lying in the 
(111) plane (see Ref. [l0| for exact definitions). Insert- 
ing Hc,3(0) = 0.05 T and the calculated values of D x x 
given in Ref. in Eqn. (f2"3"| . we obtained ( Ji + J 2 )/fc,B = 
— 5.2mK. Unlike for the other EuX compounds, we were 
not able to derive any further reliable condition that 
would allow to determine 3\ and J 2 independently of the 
high pressure data in EuSe. 

Since our Heisenberg spin model can only generate FM 
and AFM II orderings, we are only able to reproduce 
the behavior of Tc as a function of lattice constant in 
the strain-induced ferromagnetic phase of EuSe at pres- 
sures above 0.5 GPa. In this region, third nearest neigh- 
bor exchange and dipolar interactions are not expected 
to contribute significantly to the ferromagnetic ordering. 
Therefore, we calculated the ferromagnetic ordering tem- 
perature as functions of the NN and NNN Heisenberg ex- 
change interaction in the range 0.2 K < J± < 1.0 K and 
-0.175K > J 2 > -0.4K with J x > \J 2 \. As shown in 
Fig. [Tl] Tc(Ji, J 2) is slightly nonlinear in both J\ and 
J 2 , but approaches the relation given in Eqn. (|21j) for 
J 1 ^> — J 2- We also simulated the AFM II to paramag- 
netic transitions for J\ < |</ 2 | in the interval 0.13K < 
Ji < 0.18 K and -0.15 K > J 2 > -0.2 K. The calculated 
critical temperature Tjv as a function of the exchange 
constants 3\ and J 2 are shown in Figs. [TlTa) and (b) 
as open squares. As in the case of EuTc, Tjv(Ji, J 2 ) is 
strongly nonlinear and both Tq{J\, J2) and T/v(Ji, J 2 ) 
were approximated by second order polynomials as given 
in Eqn. (fl"9)) . 

Inserting the magnetic Griineisen law of Eqn. (|18p 

^i( r i) = Ji(a>/V2) and ^2 (^2) = J2 (a) into the obtained 
Tc(Ji, J 2 ) dependence, the distance dependence of the 
EuSe exchange constants Ji(a) was again obtained by 
fitting the calculated Tc(a) to the experimental Curie 
points of EuSe under hydrostatic pressure (filled sym- 
bols). Other than in the preceding sections, not only 
the Griineisen exponents ri\ and n 2 but also the ambient 
pressure exchange constants Ji ; o and J2,o were used as 
adjustable parameters in the fit routine, only restricted 
by the condition (Ji,o + J 2 ,o)/kB = — 5.2mK, obtained 
from the critical field Hc,3(0) as described above. Un- 
like the situation in EuO and EuS, the bahavior of J 2 (a) 
influences the magnetic ordering considerably in the re- 
gion close to do- Eventually, J^oA's = 0.223 ± 0.016 K, 



mate of ( J\ + J2) can be calculated using the mean field 
relation 



(23) 

I 

Ji,o/k B = -0.228 ± 0.016K, n x = 24.9 ± 1.8 and n 2 = 
12.2 ± 6.0 are obtained by the fit. As is demonstrated by 
the solid line in Fig. llOi with these parameters the exper- 
imental Tc{a) data are precisely reproduced. The depen- 
dence of Ji and J 2 versus lattice constant are depicted as 
dashed, respectively, dash-dotted lines in Fig. [TDJ At the 
bulk EuSe lattice constant of ao = 6.191 A, the calcu- 
lated J 2 is slightly larger in absolute value than J\ . This 
changes drastically as the lattice constant is reduced un- 
der hydrostatic pressure, with J\ crossing J 2 already at 
low hydrostatic strain and J\ becoming the dominant 
exchange interaction for a < 6.15 A. The Griineisen ex- 
ponents for the NN and the NNN exchange interaction 
are again in good agreement with the results obtained for 
EuTe, EuO and EuS. 

To further justify our results on EuSe, we substi- 
tuted the obtained J\(a) and J 2 (a) into the theoreti- 
cal T/v(Ji, J 2) AFM II to paramagnetic Neel function, 
obtained from the fit of Eqn. (fT^|) to the squares in 
Fig. EJ Extrapolating T N ( Ji(a), J 2 (a)) to the bulk 
lattice constant ao = 6.191 A of EuSe, we obtained 
^jvafm n( a o) = 4.5 K, which is, as expected, above the 
experimentally observed AFM II to FiM transition tem- 
perature of around 2K2*ii, but below the AFM I Neel 
point of 4.7 ± 0.1 K&il^ (open symbol in Fig. [TUJ • That 
the calculated AFM II to paramagnetic transition tem- 
perature of Tjy ^ FM n (eto) is very close to the experimen- 
tally observed AFM I to paramagnetic Ia^afm i is also 
expected in mean field theory, where 

^A A FMi(4i 2 ) = ASiS+l).!, (24) 
^Ii M „(Ji,J2) = -45(5+ 1)J 2 . (25) 

Thus, Ta^afm 1 and Ia^afm ii are nearly equal when J\ 
and J 2 are almost equal in strength. Moreover, it can 
easily be shown, that taking third-nearest neighbor or 
biquadratic exchange terms into account would not fa- 
vor cither of the two antifcrromagnctic ordering types 
(at least in the mean field approximation) and the in- 
fluence of dipolar interactions on the ordering temper- 
atures is typically of the order of less than one Kelvin 
(see, e.g., chap. 4 of Ref. [H). In addition, Lechncr et 
al.'sii results on biaxially strained EuSe indicate, that 
the AFM I ordering observed in unstrained EuSe with 
Tjv.afm 1 ~ 4.7 K is energetically only slightly lower than 
the Neel point of the AFM II ordering and small biaxial 
strain produces a transition from one to the other. 
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Nearest neighbor exchange interaction Jj (K) Next nearest neighbor exchange interaction J 2 (K) 



FIG. 11: Monte Carlo calculated critical temperatures of a bulk fee Heisenberg system as a function of the strength of (a) the 
nearest and (b) next nearest neighbor exchange interactions Ji and Ji, respectively, showing ferromagnetic to paramagnetic 
transitions if Ji > IJ2I and antiferromagnetic to paramagnetic transitions if Ji < |«7a | - The parameter range corresponds to 
EuSe under hydrostatic pressures between and around 15 GPa. Solid lines represent polynomial fit functions, which are 
quadratic in Ji and J2. 





EuTe 


EuSe 


EuS 


EuO 


ao (A) 


6.598 


6.191 


5.956 


5.144 


T N (a ) (K) 


9.85 


4.7 






T c (a ) (K) 






16.6 


69.15 


Ji,o (K) 


0.192 


0.223 


0.427 


1.169 


m 


20.6 


24.9 


22.4 


19.6 


ai (A- 1 ) 


4.56 


5.99 


5.47 


5.44 


J2,0 (K) 


-0.313 


-0.228 


-0.306 


-0.414 


n 2 


10.4 


12.2 


10*> 


10*) 


a 2 (A- 1 ) 


1.63 


2.31 


1.68*) 


1.94*) 



TABLE II: Distance dependence of nearest (subscript 1) and 
next-nearest neighbor (subscript 2) exchange interactions of 
europium chalcogenides under hydrostatic pressure as ob- 
tained in our Monte Carlo study. For hydrostatic strains 
of typically less than ~ 10%, the exchange interactions 
can equally well be described by either magnetic Griineisen 
laws, Ji — Jifi(ri/ro.i)~ ni or by simple exponential laws, 



Ji = Ji,o exp 
sumptions 



[—ai (r 



-values labelled by *' are as- 



VII. DISCUSSION 

The exchange constants of all four FiuX compounds 
and their dependence on the lattice parameter obtained 
by our MC analysis are summarized in Tab. [TT] and 
Fig. [T2] Evidently, in all cases, the bulk exchange con- 
stants Ji,o and J 2i o ar( 3 almost a factor of two larger 
compared to the values reported in previous works, which 
in their analysis did not take spin fluctuations into ac- 
count. More importantly, we have found that the mag- 
netic properties and critical phase transition tempera- 
tures of all EuX compounds as a function of hydro- 
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FIG. 12: (color online) Exchange interactions in europium 
chalcogenides as functions of the lattice parameter. 



static strain, i.e., lattice parameter, can be consistently 
described by the magnetic Griineisen power law depen- 
dence Ji(ri) = Ji t o(r.i/ri_o)~ ni with characteristic power 
law exponents of n\ w 10 for the NN exchange J\ and of 
77,2 ~ 20 for the NNN exchange J2 as indicated in Tab.HIl 
The corresponding dependence of the exchange inter- 
actions on the lattice constant J\(a) and J2(a) for all 
compounds are represented as dashed lines in Fig. 1121 
Evidently, neither Ji(a) nor J2(a) are continuous func- 
tions over different members of the EuX family. There- 
fore, the effect of the substitution of the anion elements 
cannot be simplified to a variation in the lattice constant 
alone, as was already noted in the previous work of Gon- 
charenko et al^. In fact, as shown by Fig. [T^l for a 
fixed lattice constant the absolute values of the exchange 
constants J\ and Ji are always much larger for the com- 
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pound with larger anion clement, i.e., 

|jf uTc (a)| » |jf uSo (a)| » |jf uS (a)| » |jf u °(a)| (26) 

For this reason, at a given lattice parameter the ordering 
temperature is always significantly larger in EuTe com- 
pared to EuSe, EuS and EuO. This is consistent with 
the extrapolation of the Tc(a) data obtained from ex- 
periments. 

Apart from the empirical magnetic Gruneisen law of 
Eqn. (fT8|) as the inter-ion distance dependence of the ex- 
change interactions in the EuAs, we have also explored, 
if different functional dependencies might reproduce the 
experimental ordering temperatures as functions of the 
lattice constant as well. Due to the strong superlincar 
increase of the ferromagnetic ordering temperature with 
decreasing lattice constant and the fact that Tc(Ji, J2) 
behaves very close to linear in all of the EuA compounds, 
a linear distance dependence of the exchange interactions 
can be categorically ruled out. One other empirical form 
of J(r), which could be expected from the nature of quan- 
tum mechanical two-electron multicenter integrals, e.g., 
the Hcitlcr-London approach to solve the Schrodingcr 
equation for the hydrogen molecule^!, is a simple expo- 
nential distance law of the form 

J(r) = Jo exp [-a (r - r )] (27) 
= J exp [-a r (r/r - 1)] , (28) 

where J(ro) = Jq is again the exchange constant un- 
der normal condition. Applying this exponential law 
in the simulated critical temperatures Tjv(Ji,J 2 ) and 
Tc(J 1; J 2 ) and fitting the experimental critical points of 
the EuAs as functions of the lattice constant, it turns 
out, that the results for Ji(a) and J 2 (a) are practically 
indistinguishable from those received by exploring the 
Grimciscn power law. In essence, the fitted T/v(a) and 
Tc{a) curves and the Ji(a) and J 2 (a) distance depen- 
dencies of the exchange interactions coincide almost ex- 
actly with those presented in Fig. Q21 The reason for 
this unambiguity is that the variations in the interatomic 
distances achievable by hydrostatic strain are generally 
too small (< 8%) to be able to definitely single out be- 
tween the two functional dependencies. However, while 
the distance dependencies of the exchange interaction for 
the different members of the EuA family show consis- 
tant Gruneisen exponents for the NN and the NNN 
exchange interactions, such similarities can not be ex- 
pected for the scaling factors ai in the simple exponen- 
tial description. This follows from the fact that only the 
magnetic Gruneisen law of Eqn. (|18[) is defined in terms 
of relative changes of the lattice parameter r/rg. This, 
however, can also be achieved by redefining the exponen- 
tial distance law as a function of r/r$ (Eqn. ([28)) ). Then 
a ro represents a dimensionless scaling factor that can be 
compared for the different members of the EuA family. 
The calculated scaling factors a, are given together with 
the Hi and the Jo,i in Tab. Inland are in the range of 4.5- 
6 A -1 for 011 and 1.6— 2.3 A -1 for a 2 . As it turns out, the 



values for a,i ro are approximately equal to the respective 
UiS for all EuA compounds. This becomes clear from a 
more mathematical point of view. Using the relation^ 

dJf\ 

a = -^{ , (29) 
J(r) 

and requiring that both J(r) and its derivative with re- 
spect to r have to match at ro yields the relation n = a ro. 
Thus, we have shown that for small deviations of r from 
ro the power law and the exponential dependence are 
nearly equivalent when a is chosen as a = n/r . 

The theoretical derivation of analytic scaling laws for 
the distance dependence of the exchange interactions is 
far from trivial and estimates can only be obtained as far 
as the exchange mechanisms are understood. According 
to Kasuya£ the most important contribution to the NN 
exchange interaction consists of a virtual excitation of a 
Eu 2+ 4/ electron to the 5d state of a NN cation and a 
subsequent intra-atomic d-f exchange. For this type of 
exchange, a distance dependence Ji(ri) ~ cxp(— 8ri/ro) 
is considered. The NNN exchange may consist of sev- 
eral competing components, which are considered to in- 
volve excitations of the anion p electrons to neighboring 
cation 5d states. To our knowledge Kasuya did not give 
an estimate for the distance dependence of the NNN su- 
perexchange. Lee and Liu£i on the other hand proposed 
intcrband exchange mechanisms for both J\ and J 2 where 
the exchange of the localized 4/ moments is mediated by 
virtual excitations of chalcogenide-valence band p elec- 
trons into the empty Eu 2+ 5d conduction bands, together 
with a subsequent interband exchange of the d electron 
(p hole) with the localized 4/ electrons. In this semicon- 
ductor analogue of the Ruderman-Kittcl-Kasuya-Yosida 
(RKKY) interaction J(r) ~ r~ A is considered for the 
distance dependencies of the NN and NNN exchange in- 
teractions. Both estimates for the distance dependence 
of the exchange interactions more or less support the 
picture^^^ of the exchange interactions being continu- 
ous functions across different members of the EuA fam- 
ily. Our analysis clearly shows that the exchange inter- 
actions are varying much stronger with the interatomic 
distances than previously assumed and that the exchange 
of the anion element can not be attributed solely to shifts 
in the lattice parameter. With a Gruneisen exponent of 
ri\ w 20, there is a particularly strong dependence of 
the NN exchange J\ as function of the interatomic dis- 
tances and the exponent n 2 « 10 for the distance depen- 
dence of the NNN exchange interaction J 2 coincides with 
Bloch's^i empirical law for the volume dependence of the 
superexchange in magnetic solids. 

The possibility to integrate EuO with Si and Ga N 31 i 32 
together with the fact that Tq can be increased by 
dopin g 33 ' 34 and hydrostatic strai n 44 ! 45 ' 90 to temperatures 
up to 200 K led to a renewed interest in the ferromagnetic 
europium monochalcogenides as possible materials for fu- 
ture spintronic devices. As hydrostatic pressure is not an 
option for practical applications, epitaxial strain has been 
suggested as an alternative way to increase the ferromag- 
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netic ordering temperature in EuO^I. Ingle and Elfimov 
showed in their ab initio study that biaxial compressive 
strain increases Tq in EuO similarly to the situation in 
EuSeii and EuT o 52 ' 53 , where the antiferromagnetic or- 
dering temperature is increased with a reduction in the 
in-plane lattice constant of ultrathin (111) oriented epi- 
layers. As shown by Ref . I53L by proper adjustment of the 
exchange constants, the behavior of epitaxially strained 
EuTe layers can well be described by our Monte Carlo 
calculations. 



VIII. CONCLUSION 

In summary, we have applied the Monte Carlo method 
to determine the exchange integrals in the EuA com- 
pounds. To this end, we have determined the gen- 
eral dependences of the magnetic ordering tempera- 
tures, i.e., ferromagnetic Tc(Ji, J2) and antiferromag- 
netic Tjv(Ji, J2) as functions of the NN and NNN ex- 
change interactions of a system of classical Heisenberg 
spins at the sites of an fee lattice. This was subsequently 
applied to determine the distance dependence of the NN 
and NNN exchange interactions from hydrostatic pres- 



sure experiments based on the magnetic Griineisen law- 5 - 7 - 
Ji{ r i) ~ r i~™% where r, denotes the interatomic distances 
between neighboring Eu ions. It turns out that the dis- 
tance dependences of the exchange interactions J\{r\) 
and 72(72) of the different members of the EuA fam- 
ily can be consistently described by Griineisen exponents 
m ~ 20 and n-i ~ 10, where the latter conforms with 
Bloch's empirical 10/3 law for the volume dependence 
of the superexchange interaction 5 - 7 -. The strong depen- 
dence of the exchange constants on the lattice parameter 
provides room for substantially increasing the magnetic 
ordering temperatures in strained heteroepitaxial struc- 
tures, which is an important prerequisite for device ap- 
plications. 
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